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Abstract 

The mechanism of cohectivity coexisting with chaos in a finite system of 
strongly interacting fermions is investigated. The complex spectra are rep- 
resented in the basis of two-particle two-hole states describing the nuclear 
double-charge exchange modes in ^^Ca. An example of J'^ = 0~ excitations 
shows that the residual interaction, which generically implies chaotic behav- 
ior, under certain specific and well identified conditions may create strong 
transitions, even much stronger than those corresponding to a pure mean- 
field picture. Such an effect results from correlations among the off-diagonal 
matrix elements, is connected with locally reduced density of states and a 
local minimum in the information entropy. 

PACS numbers: 05.45.+b, 05.30.Fk, 21.60.Ev, 24.30. Cz 
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I. INTRODUCTION 



The concept of the random matrix theory (RMT) proves very fruitful in approaching 
complex quantum systems and in addressing the question of how classical chaos manifests 
itself on the quantum level. Chaos is essentially a generic property of complex systems 
such as atomic nuclei 0, many electron atoms 0, molecules |^ or disordered mesoscopic 
systems |^ and this finds evidence in a broad applicability of RMT to describe level fluc- 
tuations 1^. Even many aspects of quantum chromodynamics are consistent with chiral 
RMT [0]. Similarly, however, as in most physically interesting cases where classical chaos 
is not just a hard billiard-type chaos, the pure RMT cannot account for the full richness of 
quantum phenomena connected with complexity. As an example one can mention the sign 
correlations [§] for parity nonconserving effects in compound nuclei, even though it was 
the physics of compound nuclei which lead Wigner [|lOl to postulate the Gaussian orthogonal 
ensemble (GOE) of random matrices as an appropriate global frame. Explicit microscopic 
approaches in terms of the full shell model diagonalization, either in nuclear [0,|1^ or atomic 
physics [0], show perfect agreement with GOE when looking at the local level fiuctations 
measured in terms of the nearest neighbor spacing distribution and the A3 statistics, but 
significant deviations take place on the level of wave functions. This originates from the 
two-body nature of interaction which reduces the number of independent parameters and 
preserves certain correlations among the matrix elements. In order to account for this type 



of correlations a two-body random interaction model has been introduced ||I2| and its sta- 



tistical properties investigated in detail ||T3|. Still however, such models may not properly 
account for correlations which originate from geometry of a problem and which, in some 
cases, may turn out significant. 

Another characteristics connected with complexity, even more interesting and important 
from the practical point of view, is collectivity. It means a cooperation, and thus the 
coupling between the different degrees of freedom in order to generate a coherent signal 
in response to an external perturbation. Consequently, even though the real collectivity 
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implies a highly ordered behavior it involves effects beyond the mean-field - the most regular 
part of the many-body Hamiltonian. At the same time the effects beyond the mean 
field are responsible for the GOE fiuctuation properties. Therefore, in a sense, these two 
seemingly contradictory phenomena, chaos and collectivity, may have to go in parallel. Also 
on classical level collectivity is a nonlinear cooperative effect which results from the coupling 
between different degrees of freedom. 

In general, the shell model type approaches are based on diagonalization of the full 
many-body Hamiltonian in the basis spanned by all possible n-particle - n-hole (npnh) 
configurations generated by the mean-field. For practical reasons, especially when large 
energy intervals are involved, as for instance in the case of nuclear giant resonances, one 
truncates this hierarchy of configurations up to 2p2h [|^. Interestingly, due to a sufficiently 
large density of states relative to the strength of the residual interaction |jl6[, local level 



ffuctuations characteristic of GOE appear [|T^ to take place for the nuclear Hamiltonian 



acting already in the space of 2p2h states and this is a crucial element for an appropriate 
description of the giant resonace decay properties [|TB|. The giant resonaces are however 
excited by one-body operators which directly probe the Iplh components of the nuclear wave 
function. The 2p2h states only form the background which determines a decay-law. There 
exist, however, very interesting physical processes, represented by two-body (two-phonon) 
external operators, which directly couple the ground state to the space of 2p2h states. In 
view of the above mentioned local GOE ffuctuations giving evidence for a significant amount 
of chaotic dynamics already in the 2p2h space, the question of a possible coherent response or 
collectivity under such conditions is a very intriguing one and of interest for many branches 
of physics. 

II. MODEL 

We start from the Hamiltonian which, in second quantized form, reads as 
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H = Y1 ^i^i^i i ^ Vij^kiajo'laiak- (1) 

i ij,kl 

The first term denotes the mean field while the second term is the residual interaction with 
antisymmetrized matrix elements Vij^ki- Diagonalizing this Hamiltonian in the subspace of 
2p2h states 

\2) = alyp^tth^ahAO) (2) 

yields the eigenenergies En and the corresponding eigenvectors |n) = S2C2I2). For realis- 
tic nuclear interactions the spectral fiuctuations of {En} typically coincide with those of 
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The general form of matrix elements for the two-body residual interaction V between 
2p2h-states is given by 

^piP2hih2,p[p'2h[h'2 

= Sp^p'^5p2p'^Vh[h!^hih2 + Shih[Sh2h'2'^^piP2p[p'2 (3) 

+ a{Pi,P2)a{hi, h2)a{p[,p2)a{h[, h2)6p^p'^6h2hiVpih[hip[, 

where a{i,j) denotes the antisymmetrizer between i and j. The consecutive terms in this 
expression are responsible for hole-hole, particle-particle and particle-hole interactions, re- 
spectively, while the remaining pair of states in each case are spectators represented by the 
6ij functions. These functions set a significant fraction of the matrix elements to zero which 
may lead to correlations. Fig. 1 illustrates the corresponding structure in diagrammatical 
representation. Further correlations may originate from the fact that many nonzero matrix 
elements relate to each other only by the geometrical factors due to the angular momentum 
coupling algebra. 

In response to an external field Ea a state 

\E^) = F^\0)=J2{n\F^\0)\n) (4) 

n 

is excited. The two-phonon operator E^ can be represented as 



Fa = [U ® flU (5) 

where fp and /y denote the single-phonon operators whose quantum numbers (5 and 7 are 
coupled to form a. The state \Fa) determines the strength function 

Sf^{E) = Y.SfMKE-E^), (6) 

n 

where 

SpM = MFam\ (7) 

In the unperturbed basis of states |2) the transition strength Sp^in) to the state |n) can be 
expressed as 

= E Ic2ri(2|4|0)|^ + E «(0|^i|2')(2|4|0) 

2 2^2' 

^S'pAn) + S]i{n). (8) 

The second equahty defines the diagonal {Sp^n)) and off-diagonal {Sj^{n)) contributions to 
the transition strength at energy E^- The second component includes many more terms and 
it is this component which potentially is able to induce collectivity, i.e. a strong transition 
to energy E^- Two elements are however required: (i) a state \n) must involve sufficiently 
many expansion coefficients C2 over the unperturbed states |2) which carry the strength 
((2|Fa|0) 7^ 0) and this is equivalent to at least local mixing, but at the same time (ii) sign 
correlations among these expansion coefficients should take place so that the different terms 
do not cancel out. 

Optimal circumstances for the second condition to be fulfilled read: 

c^~(0|F«|2). (9) 

This may occur if the interaction matrix elements can be represented by a sum of separable 
terms of the multipole-multipole type 



M 

v.,,M = E QWm (10) 



with Q\- ~ The success of the Brown-Bolsterh schematic model |]T^ in indicating 

the mechanism of collectivity on the Iplh level points to an approximate validity of such a 
representation and its formal justification comes from the multipole expansion of the residual 
interaction. The structure of the Hamiltonian matrix in the Iplh-subspace is then usually 
dominated by few multipoles. Collectivity can then be viewed as an edge effect connected 
with appearance of a dominating component in the Hamiltonian matrix and the rank M of 
this component is significantly lower (unity in case of the Brown-Bolsterli model) than the 
size of the matrix. This rank specifies a number of the prevailing states whose expansion 
coefficients predominantly are functions of Q'^. In general, on the 2p2h level a multipole 
structure of the interaction enters the corresponding matrix elements in a more complicated 
way. However, due to two-body nature of the nuclear interaction which reduces its 2p2h 
matrix elements to combinations of the ones representing the particle-particle, hole-hole and 
particle-hole interactions [|l5l the separability may become effective also on the 2p2h level 
although conditions are expected to be more restrictive. On the other hand the 2p2h space 
offers many more unperturbed transitions to form a collective state and the net effect may 
still appear significant. 

For quantitative discussion presented below we choose the ^^Ca nucleus, specify the mean 
field part of the Hamiltonian (|l|) in terms of a local Woods-Saxon potential including the 



Coulomb interaction and adopt the density- dependent zero-range interaction of ref. ||2^ as a 
residual interaction (after correcting for a misprint in the density functional: Rq = 1.16A^^^). 
Since we want to inspect the higher energy region at least three mean field shells on both 
sides of the Fermi surface have to be used to generate the unperturbed 2p2h states as a 
basis for diagonalization of the full Hamiltonian (|I|). Typically, the number of such states 
is very large and this kind of calculation can be kept under full numerical control only for 
selected excitations of the lowest multipolarity. Among various nuclear excitation modes 
which can be considered in this context the double-charge exchange (DCX) processes are 
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of special interest. These modes, excited in (vr"*",??") reactions involve at least two 
nucleons within the nucleus and give rise to a sharp peak at around 50 MeV in the forward 
cross section. They are thus located in the energy region of the high density of 2p2h states 
which points to the importance of coherence effects among those states. Consequently, the 
present investigation may also appear helpful in studying the mechanism of DCX reactions 



and in separating the suggested dibaryon contribution from the conventional effects [[23 
For all these reasons we perform a systematic study of the DCX J'^ = 0~ states. Our model 
space then develops N=2286 2p2h states. There are still several possibilities of exciting such 
a double-phonon mode represented by the operator Fa out of the two single-phonons f/3 and 
of opposite parity. For definitness we choose 

fp = rY,T^ (11) 

and 

/^ = r2[F2®a]i+r_. (12) 

The first of these operators corresponds to the ITiuo dipole and the second to 2huj spin- 
quadrupole excitation. The resulting two-phonon mode thus operates on a level of 3^a; 
excitations. Formulas needed to express the angular momentum coupled form of the above 



one- and two-body operators can be found, for instance, in the appendix of ref. |24 



III. RESULTS AND DISCUSSION 

The results of calculations are presented in Fig. 2. As one can see, including the residual 
interaction (part (b)) induces a strong transition at 51.1 MeV. This transition is stronger 
by almost a factor of 2 than any of the unperturbed (part (a)) transitions even though it is 
shifted to a significantly higher (~ 10 MeV) energy. This is also a very collective transition. 
About 95% of the corresponding strength originates from 5'|^^(r^,), as comparison between 
parts (b) and (c) of Fig. 2 indicates. This whole effect is due to particle-hole type matrix 
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elements ((c) in Fig. 1). Discarding diagrams (a) and (b) produces no significant difference. 



The degree of mixing can be quantified, for instance, in terms of tlie information entropy |25 



I {n) = -T^iPi In Pi] Pi = |c"| (13) 

of an eigenvector \n) in tlie basis (part (d)). Interestingly, the system finds preferential 
conditions for creating the most collective state in the energy region of local minimum in 
/(n). Our following discussion is supposed to shead more light on this issue. 

As shown in Fig. 3(a) our Hamiltonian matrix displays a band-like structure with spots of 
the significant matrix elements inside. This together with a nonuniform energy distribution 
Pu{E) of the unperturbed 2p2h states (Fig. 3(b)), which is a trace of the shell structure of the 
single particle states, characteristic of many other mesoscopic systems [^], sizably suppresses 



the range of mixing and locally supports conditions for the edge effect to occur in the energy 
region of the minimum in p{E). A comparison with Fig. 2(b) shows that the collective state 
is located at about this region. Moreover, the minimum survives diagonalization {pp{E) in 
Fig. 2(c)) and all the above features are consistent with the effective band range 



{AE^r = E(^^^ - HnYny Hi (14) 

shown in Fig. 2(d). 

Further quantification of the character of mixing between the unperturbed states is 
documented in Fig. 4. The distribution P{H) of off-diagonal matrix elements (a) is not 
Gaussian but of the following type: 

P{H) = a\H\^exp{-\H\/c) (15) 

This indicates the presence of the dominating multipole-multipole components in the in- 
teraction [|2|,|^. An interesting feature is the asymmetry between the positive and negative 
valued matrix elements (see parameters in the caption to Fig. 4(a)). The positive matrix ele- 
ments are more abundant which expresses further correlations among them and the fact that 
the interaction is predominantly repulsive for the mode considered. Significant reduction of 
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dimensionality is also indicated by the distribution of eigenvalues of the residual interaction. 
As shown in part (b) of Fig. 4 the majority of these eigenvalues is concentrated around 
zero and thus constitute approximate zero modes of that part of the Hamiltonian. We also 
would like to note at this point, without showing the results explicitly, that similar analysis 
on the Iplh level using appropriately larger model spaces (for better statistics) shows an 
even larger fraction of such zero modes. This is due to the fact that in the Iplh-space the 
multipole-multipole structure of the interaction manifests itself in a more transparent way. 

Appearance of a strong transition at certain energy En means that the structure of the 
Hamiltonian matrix of the residual interaction, at least locally at around that particular 



energy, is governed by a component of the type as specified by eq. (|TOD with a small number 
of terms (M <^ A^) including, of course, the ones which coincide with an external field. In 
a pure case a structure like this causes an energy gap between the collective state and the 
remaining states. In the present case, of its only local nature, one expects a local minimum 
in the density of states in the vicinity of the collective state. Indeed, as can be seen by a 
careful inspection of Fig. 2(b) versus Fig. 3(c) any stronger transition is located in such a 
minimum whose range typically extends over an energy interval of the order of 0.5 MeV. 
Even relatively weak transitions are asigned their own minima. Moreover, as we have verified 
in certain selected cases, many other minima in the density of states which are not occupied 
by the above specified transitions turn out to be filled in by the DCX J'^ = 0~ transitions 
connected with other combinations of two one-phonon operators (for instance = r'^Y2T- 
and fj = r[Yi cr]2-r_). 

A reduction of the rank (real dimensionality) of the Hamiltonian matrix evidenced above 
is also consistent with the observed minimum in the information entropy (Fig. 2(d)). Sim- 
ply, in the relevant energy region there are fewer free parameters and this sets additional 
constraints on the degree of mixing and thus on the amount of chaos. As a chaos related 
characteristics we take the spectral rigidity measured in terms of the A3 statistics We 
find this measure more appropriate for studying various local subtleties of mixing than the 
nearest neighbor spacing (NNS) distribution because for a smaller number of states the latter 
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sooner becomes contaminated by strong fluctuations. Indeed, the spectral rigidity (Fig. 5) 
detects differences in the level repulsion inside the string of eigenvalues (il) covering the flrst 
maximum in ppi^u){E) (35.2-44.5 MeV, 400 states starting from n = 351 up to n = 750) and 
the one (12) covering the minimum and thus including the collective state (44.5-52.1 MeV, 
400 states from n = 751 to n = 1150). The deviation from GOE is more signiflcant in 12 
which, similarly as I{n), signals a more regular dynamics in the vicinity of the collective 
state {n = 1089). 

Conditions corresponding to the actual Hamiltonian are not the most optimal ones from 
the point of view of the collectivity of our J'^ = 0^ DCX excitation. By multiplying the 
residual interaction by a factor of (? = 0.7 we obtain a picture as shown in Fig. 6. Now the 
transition located at 48.6 MeV is another factor of 2 stronger than before and, again, all 
signiflcant transitions are situated in the local minima of p{E) and in the overall minimum of 
the information entropy. A too severe decrease of g will eventually bring all the transitions 
to their mean-fleld values. This transition is however not just linear. Here we seem to 
be facing a competition of the two elements. One is the residual interaction which must be 
sufficiently strong to correlate many states but the other one is condition for the eddge effect 
to occur. As a result, even for g = 0.35 we still obtain very strong transitions, apparently 
due to the fact that the interaction strength is such that the unperturbed transitions are 
moved just to the absolute minimum in p{E). 

The range of values of a multiplication factor which produces this kind of picture is rather 
narrow and this feature of collectivity resembles a classical phenomenon of the stochastic res- 



onance 



p8| . It is relatively easy to completely destroy so strong transitions. By multiplying 
the residual interaction by a factor of g = 2.5 (which is equivalent to increasing the density 
of states) the strength distribution displays a form as shown in the lowest panel of Fig. 6. 
This strength remains largely localized in energy but the distribution of the correspond- 
ing Spin) (Fig. 6) does not deviate much from the Porter-Thomas (P-T) distribution ||2^ 



P{s) = (27rs)-^/2gxp(-s/2) characteristic of GOE, even though correlations among the 
matrix elements are the same as before. Interestingly, even here the larger transitions are 
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located in their own small minima in p{E). Further increase of the multiplication factor 
may again produce some transitions which are more collective than those allowed by P-T. 
In particular, starting from values ~4 some new strong collective transitions appear at the 
upper edge of the whole spectrum. 

To illustrate statistics of the transition strength versus P-T distribution we use a measure 
introduced in ref. Consequently, for all the cases considered above we calculate the total 
number N of transitions of magnitude smaller than a given threshold value Sth, as a function 
of Sth- Since the number of large components relative to the small ones is of primary interest 
in the present study and in order to set the same scale when comparing different cases we, 
in addition, in each case independently, devide all the transitions S{n) by the corresponding 
maximum value of S{n). After that Smax{i^) = 1 in each case. Consistently, the RMT 
limit of this measure is then drawn from the cummulative P-T distribution and this limit is 
indicated by the solid line in Fig. 7. In the log-log scale this limit develops a long straight 
line segment with the slope of 0.5 which reflects the dominant role of the preexponential 
factor in P-T distribution at smaller transitions. As it can be seen from this Figure, 

the g = 2.5 case is very close to this limit. But, interestingly, even g = 0.35 tends to the 
same slope when probing the region of small transitions which means that such transitions 
are consistent with GOE. Only the unperturbed case is distinct in this sense. 

Finally, Fig. 8 shows the transition strength distribution of the 'constituent' single- 
phonon modes specified by eqs.(0) and (p!2D, respectively, in their own Iplh sectors and the 
same model space of single-particle states is used. There are 28 (single) charge exchange 
J'^ = 1~ and 25 J'" = Iplh states in this space. As before, the ^f-factors reflect the 
strength of the residual interaction relative to the original one. The results collected in this 
Figure provide further evidence that collectivity observed on the 2p2h level is not acciden- 
tal. It can always be traced back to collectivity of the corresponding single-phonon modes in 
their own subspaces. Consistently with our previous discussion, this correspondence cannot 
however be expressed simply in terms of proportionality. For instance, the single-phonon 
transitions for g = 0.35 are significantly weaker than for g = 1 while the opposite applies 
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to the resulting two-phonon mode. Recent study of ref. |3^, even though based on a much 
simpler model, also shows that characteristics of the two-phonon mode (double giant dipole) 
are much more sensitive to detailed form of the Hamiltonian than those of the corresponding 
single-phonon modes. 

Taken together, a real collectivity, by which we mean a transition stronger than those 
generated by the mean field, is a very subtle effect and is not a generic property of the 
complex spectra. Its appearance, as it happens for one of the components of the J'^ = 
0^ DCX excitations considered here, involves several elements like correlations among the 
matrix elements, nonuniformities in the distribution of states and a proper matching of the 
interaction strength to an initial (unperturbed) location of the transition strength relative 
to the scale of nonuniformities in the distribution of states. If present, a collective state 
is then located in the region of more regular dynamics characterized by lower information 
entropy, more sizable deviations from GOE of the level fluctuations and local minima in the 
density of states. This later effect can thus potentially be used in experimental studies as an 
extra criterion to detect collectivity. We also would like to point out that these aspects of 
collectivity parallel an analogous property hypothesised for living organisms [^] and stating 
that collectivity is a phenomenon occuring at the border between chaos and regularity. 

This work was supported in part by Polish KBN Grant No. 2 P03B 140 10 and by the 
German-Polish scientific exchange program. 
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FIGURE CAPTIONS 



Fig. 1. Diagrammatic representation of the two-body matrix elements in the space 2p2h 
states with exphcit indication of the angular momentum couphng scheme. The consequtive 
terms represent hole-hole, particle-particle and particle-hole interactions, respectively. 
Fig. 2. (a) The unperturbed transition-strength distribution in ^^Ca for the J'^ = 0~ DCX 
excitation involving the single-phonon dipole and 2huj spin-quadrupole modes, (b) The same 
as (a) but after including the residual interaction, (c) Sp^ components of the transition- 
strength as defined by eq. (§) (notice different scale), (d) The information entropy of the 
states \n) in the unperturbed basis. The dashed line indicates the GOE limit (ln(0.48A^)). 
Fig. 3. (a) Structure of the Hamiltonian matrix for the J'" = 0^ DCX states. The states 
are here labeled by energies, ordered in ascending order and the matrix elements Hn. > 0.1 
are indicated by the dots, (b) density of the unperturbed 2p2h-states. (c) density of states 
after the diagonalization. (d) energy range of interaction between the unperturbed states. 
Fig. 4. (a) Distribution of off-diagonal matrix elements between the J'^ = 0~ DCX states 
(histogram). The solid lines indicate fit in terms of eq. 113] with the resulting parameters: 
a = 676, b = -1.21, c = 0.69 (left) and a = 692, b = -1.22, c = 0.81 (right), (b) Density of 
states corresponding to the residual interaction part of the Hamitonian (|lD. 
Fig. 5. Spectral rigidity A3(L) for eigenvalues from the two 400 states long intervals: 
n = 351 — 750 (il) and n = 751 — 1150 (12). The long-dashed line corresponds to Poisson 
level distribution and the short-dashed line to GOE. 

Fig. 6. Transition-strength distribution, density of states and the information entropy for 
the same excitation as in Fig. 1 but the residual interaction is now multiplied by a factor of 
g = 0.35 (upper part), g = 0.7 (middle part) and g = 2.5 (lower part), respectively. 
Fig. 7 The total number of transitions of given strength (properly rescaled, see text) 
below a threshold value Sth- The open crosses refer to the unperturbed case {g = 0), thick 
dots to g = 0.35, open squares to g = 0.7, filled squares to g = 1 and open triangles to 
g = 2.5. The solid line represents the same quantity determined from a Porter- Thomas 
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distribution. 

Fig. 8 Transition-strength distribution corresponding to single-phonon operators specified 
by eq. ( pTD (left) and by eq. (p^ (riglit), calculated in the space of Iplh states for various 
residual interaction multiplication factors g. 
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